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Abstract 

We discuss the energy flow of the classical gluon fields created in collisions of heavy nuclei at collider energies. 
We show how the Yang-Mills analoga of Faraday's Law and Gauss' Law predict the initial gluon flux tubes 
to expand or bend. The resulting transverse and longitudinal structure of the Poynting vector field has a rich 
phenomenology. Besides the well known radial and elliptic flow in transverse direction, classical quantum 
chromodynamics predicts a rapidity-odd transverse flow that tilts the fireball for non-central collisions, and 
it implies a characteristic flow pattern for collisions of non-symmetric systems A + B. The rapidity-odd 
transverse flow translates into a directed particle flow v% which has been observed at RHIC and LHC. The 
global flow fields in heavy ion collisions could be a powerful check for the validity of classical Yang-Mill 
dynamics in high energy collisions. 
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1. Introduction 

Hadrons and nuclei colliding at very large ener- 
gies can be used to study a novel regime of quan- 
tum chromodynamics (QCD) called the color glass 
condensate (CGC) It is characterized by a sat- 
uration of the transverse densities of gluons ~ 1/Q« 
in the initial wave functions, setting an energy scale 
Q s , the saturation scale. High occupation numbers 
permit a quasi-classical description of the gluon 
fields in the nuclei before collisions iS and at 
very early time after the collision [4|. Data from 
the Relativistic Heavy Ion Collider (RHIC) and the 
Large Hadron Collider (LHC) have given us hints 
that we have reached an energy regime where CGC 
could be the correct description of the relevant de- 
grees of freedom of early stages of nuclear collisions. 

There is a large body of evidence that a ther- 
malizcd quark gluon plasma (QGP) is eventually 
created within the first fm/c after the collision of 
nuclei at RHIC or LHC energies 0, Q. CGC as 
an effective theory can be used to describe the very 
first step in this process, from the original nuclear 
wave functions to a space-time volume filled with 



coherent and far off-equilibrium gluon fields, called 
a glasma, at a time tq ~ 1/Q S ~ 0.1 . . . 0.2 fm/c af- 
ter the collision. The dominant feature of the fields 
immediately after the collision are flux tubes of 
longitudinal chromo-electric and chromo-magnetic 
fields @, H Hi- After this initial period the fields 
have to decohere and the emer ging particles have 
to approach equlibrium 1(J 11, 12\. After equilib- 
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rium is reached at a time Tth > tq hydrodynamics 
can be used to describe t he g lobal collective mo- 
tion of the system [l3l [Til flBT 16|. The success of 
hydrodynamics in heavy ion collisions is based on 
its ability to capture key features of the collective 
flow of the thermalized QGP, foremost the radial 
transverse flow v r and the quadrupole asymmetry 
V2, called elliptic flow. 

However, flow is not a privilege of thermalized 
matter with pressure gradients. In fact the effective 
"transverse pressure" pt = T 11 = T 22 in a purely 
longitudinal classical gauge field is equal to the en- 
ergy density T 00 , and thus gradients in pressure 
would be larger by a factor 3 than in a relativis- 
tic free gas with the same energy density 0, E3]- 
Thus the amount of initial collective motion, of- 
ten termed pre-equilibrium flow, at the starting 
time T t h of the hydrodynamic evolution is not ex- 
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pected to be small. In fact various efforts have 
been made to estimate the magnitude of transverse 
flow around midrapidity from color glass or other 
dynamics 17. Hi}, although many hydrodynamics 
simulations still neglect pre-equilibrium flow. 

In this Letter we expand existing calculations 
by analyzing the full 3+1-dimensional structure 
of the initial gluon energy flow. To that end we 
will present an analytic calculation of the chromo- 
electric and chromo-magnetic fields E and B im- 
mediately after the collision of nuclei in a strictly 
classical implementation of CGC, the McLerran- 
Venugopalan (MV) model @, 0]. We then com- 
pute the energy momentum tensor of the gluon 
field as a function of space-time coordinates for 
small proper times t = \/t 2 — z 2 after the colli- 
sion. We analyze the resulting global map of en- 
ergy and momentum flow, in particular the trans- 
verse Poynting vector S l — T 0t , i = 1,2. We will 
argue that it can be decomposed into two parts, 
S l + and S' l _. S\_ is even in the space-time rapid- 
ity r] = 1/2 In [(t + z)/(t — z)] and driven by the 
gradient of the initial transverse pressure pr of the 
gluon field; in short this term behaves "hydro-like" 
and has in fact been discussed before analytically 
and numerically 0, H3, E3 ■ The novel 77-odd term 
51 arises from asymmetries between the two nu- 
clei, e.g. for collisions A+B of different nuclei or for 
symmetric A+A collision with finite impact param- 
eter b. Since it vanishes at mid-rapidity it has not 
been discussed in situations where one expects flow 
to be (almost) boost-invariant. One can show that 
despite the superficial appearance the 77-odd flow 
term Sl_ does not break the boost-invariance which 
is intrinsic to the MV model. We will argue, using 
the abelian case, that both rapidity-odd and -even 
terms are a natural consequence of Faraday's Law 
and Gauss' Law. We will conclude this Letter by 
discussing the phenomenological consequences for 
several colliding systems and impact parameters. 

2. Gluon Fields in the MV Model 

In the CGC the large momentum components 
of the nuclear wave function are assumed to be 
strongly Lorentz-contracted in the lab frame and 
can be approximated by an infinitely thin sheet. 
The large momentum quarks and gluons form an 
effective color SU(3) charge density Paif), a — 
1, . . . N 2 — 1, on that sheet, where r is a transverse 
position. They act as sources for the soft gluon 
modes which are described by the classical field A^ . 



In the McLerran-Venugopalan model the source p a 
moving along the +- light cone leads to a SU(3)- 
current = 5 t * + 5(x~)pat- in light cone gauge. 
The soft gluon modes are determined by the Yang- 
Mills equation D^F^ = J". The charge density 
Paif) varies on time scales that are large compared 
to typical time scales of the collision. It vanishes on 
average, (p) = 0, due to local color neutrality, but 
has a non-vanishing variance p which in the MV 
model is realized by a Gaussian distribution with 
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{Pa{ri)pb{r 2 )) = N2 _ ^ Sabfi(ri)5 (n - r 2 ) ■ (1) 

Thus p is the average, color- summed, squared den- 
sity of charges in the transverse plane (which differs 
from the definition in [Tsj by a factor N 2 — 1). 

The most fundamental result of the MV model is 
an expression for the gluon distribution function of 
a nucleus in light cone gauge (i,j = 1, 2) 

(Aimm = in g ( 2) 

where the UV limit (both fields evaluated at the 
same transverse position r) has been taken. Here Q 
is a UV cutoff related to the resolution with which 
the gluon distribution is probed and m is an in- 
frared scale that can be introduced as an effective 
gluon mass. This result was first calculated in Ref. 
20] for infinitely large, homogeneous nuclei, i.e. for 
constant source area density p. Of course this is 
not a realistic approximation for nuclei. We have 
checked that if variations of p(f) are permitted but 
are small on the IR length scale ~ 1/m, i.e. 



to 1 \Vp(r)\ < p(r) 



(3) 



many well-known results of the MV model, in par- 
ticular Eq. @ hold and corrections can be classified 
in an expansion in gradients m~ ly \7 l p(r). In sim- 
ple terms this indicates that the physics of color 
glass and the less known QCD dynamics at larger 
distances can be separated in a meaningful way as 
long as condition ([3]) is satisfied. We will report 
details of this calculation elsewhere 2l[. 

Given the purely transverse gauge fields A\ = 
A t a l (r)t± and A\ = A* i2 (f)t a of the nuclei before 
the collision, in their respective light cone gauge, 
one can write the longitudinal chromo-electric and 
-magnetic fields after the collision for small times r 
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as [811211,122 

E 3 
B 3 



1 + —D 2 )e + O{t^ 
1 • Oir 1 ) 



(4) 
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where D 2 = D l D t is the square of the covariant 
derivative with respect to gauge field A\ + A\ and 



Eq = ig 



A 3 A J 

-"11 ^2 



A{,A 



arc the Ion- 



Bo = ige jk 

gitudinal fields at r = 0. Since the transverse fields 
vanish at r = this has given rise to the notion 
of flux tubes of longitudinal electric and magnetic 
fields [9( . In terms of those initial longitudinal fields 
the transverse fields can be expressed as || [U [22| 

E i = ~ (sinh r]D l E + cosh 77 ( lJ D J B Q ) (6) 
B l = - T - (sinh r/D l B a - coshrj e ij D j E ) (7) 

plus terms of order r 3 and higher. 

It is straightforward to write down the initial en- 
ergy density e = T qo (t = 0) = (E%+B$)/2 and the 
initial transverse energy flux given by the Poynting 
vector T°\ The leading behavior at small r has 
four terms, two even in rj, 



? + = - cosh?? (E D l E + B D l B ) 



and two odd in rj, 



St = - sinh 77 e^' (E D 3 B q - B D 3 E ) 



(8) 



(9) 



and T 0i = S q + + S l _ . The longitudinal energy flow 
is suppressed and starts at order 0(t 2 ), 

T 03 = [ sinh2r7 ((DE f + (DB ) 2 ) 

+ 2cosh2rje ll (D i E )(D 1 B )] . (10) 

Note that the MV model is boost-invariant by def- 
inition and the ^-dependence of the energy flow 
does not violate boost-invariance. One can verify 
that the full energy momentum tensor T' 1 " satisfies 
A w '(y)T MV ,(T ? )A^'(zj) = T^{ri + y) where A(y) is 
the Lorenz boost tensor with rapidity y along the 
z-axis. In fact the rapidity-odd field is a natural 
consequence of the field equations as the next sec- 
tion will show. 




Figure 1: Two observers at z = zq and z = — zrj test 
Ampere's and Faraday's Laws with areas a 2 in the trans- 
verse plane and Gauss' Law with a cube of volume a 3 . The 
transverse fields from Ampere's and Faraday's Laws (black 
solid arrows) are the same in both cases, while the transverse 
fields from Gauss' Law (black dashed arrows) are observed 
with opposite signs. Initial longitudinal fields are indicated 
by solid grey arrows, thickness reflects field strength. 



3. An Electrodynamic Analogon 

Let us consider the following equivalent bound- 
ary value problem in classical electrodynamics. In 
the forward light cone r > we have the Maxwell 
Equations d fl F flL ' — without sources. On the light 
cone r = we demand the boundary conditions 
E(t = 0,r) = E (r)e z , B(t = 0,r) = B {r)e z , 
i.e. the initial fields are purely longitudinal. We 
also assume that those fields are related through 
transverse fields A\ and A\ as Eq — S 1 ^ A\A\ and 
Bq = e l 3 A\A 3 2 . The abelian problem for fixed initial 
conditions can in principle be solved analytically 
[H Hl| j but it will suffice here to give the solution 
order by order in powers of r as we did in the case of 
QCD. From the QCD solutions we can immediately 
conclude that the longitudinal fields in the abelian 
case are 



E' = { I + | E l} 

t 2 -z 2 



B = \ 1 



V 2 Bo , 



while the transverse fields are 

E i = ^V l E + y^Bo 
B* = ^\7*B - ^V j E 



(11) 
(12) 

(13) 
(14) 
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Figure 2: Transverse electric fields (left panels) and magnetic 
fields (right panels) at r) = (upper panels) and r\ = 1 (lower 
panels) in an abelian example for a random distribution of 
fields A^, A l 2 . The initial longitudinal fields So (left panels) 
and Eo (right panels) are indicated through the density of 
the background (lighter color = larger values). At r) = the 
fields are divergence-free and clearly following Ampere's and 
Faraday's Laws, respectively. 

for small times r, i.e. t 2 « z 2 . The cartesian coordi- 
nates permit simple checks of these solutions with 
Gauss', Ampere's and Faraday's Laws. 

There is a straight-forward interpretation of some 
aspects of these results. Let us choose, just as an 
example, a transverse position where Eq,Bq > 
and V 2 Eq, V 2 Bq < so that the longitudinal fields 
decrease away from the light cone t 2 — z 2 . Two ob- 
servers at fixed points z = zq > and z = — zq 
would observe the same electric (magnetic) flux 
through a small transverse area a 2 with an initial 
value Eqcl 2 (Boa 2 ) at t = zq which then dimin- 
ishes at the same rate V 2 E a 2 t/2 (V 2 B a 2 t/2) for 
both. Due to Ampere's (Faraday's) Law this reduc- 
tion induces magnetic (electric) fields curling with 
a negative (positive) chirality around the longitudi- 
nal fields, respectively, see Fig. [TJ 

On the other hand the same two observers at 
fixed points zq and — zq can at time t — zq count 
the electric or magnetic flux through small cubes of 
volume a 3 whose sides are aligned with the coordi- 
nate axes. One side is held at z — ±zq, while the 
opposite side is at z = ±zq qp a for the observer at 
zq or — zq respectively. In the former case the total 




Figure 3: Example for transverse flow of energy for r) = 
(left panel) and 77 = f (right panel) in the abelian example 
for the same random distribution of fields A\, A\ as in Fig. [2] 
The initial energy density T 00 is shown through the density 
of the background (lighter color = larger values). At r] = 
the flow follows the gradient in the energy density in a 
hydro-like way while away from mid-rapidity energy flow gets 
quenched in some directions and amplified in others. 

flux out of the box due to the longitudinal field is 
— Zoa 3 W 2 Eq/2 > while for the observer at — zq the 
net flux of longitudinal field has the opposite sign. 
Thus at zq a net flux of transverse field has to enter 
into the box while at — zq the same amount has to 
flow out of the box to satisfy Gauss' Law. 

To summarize, the transverse fields naturally 
have a part due to Gauss' Law with vanishing cir- 
culation (e* J 'V J ') 5 which is odd in 77, and they have a 
part due to Ampere's and Faraday's Law (and with 
different signs between the magnetic and electric 
part due to the Lenz rule) with vanishing trans- 
verse divergence (V 1 ), which is even in 77. Fig. [5] 
shows the transverse electric and magnetic fields 
for two rapidities r\ for random fields A\ and A\ 
in a sector of the transverse plane. One can check 
that this statement about transverse fields trans- 
lates directly into a matching statement about the 
transverse flow of energy since the initial transverse 
Poynting vector is linear in the transverse fields, 
T° l = ^{WBq - B j E ). Thus we have the four 
contributions already discussed in the case of QCD, 
two of them odd in 77. Fig. [3] shows the flux of en- 
ergy in the transverse plane for two rapidities for 
the same random configuration of fields A\, A\. 

4. Averaging Over Field Configurations 

Eqs. ijHJ), (H]) can be the basis for a Monte Carlo 
modeling of the energy flow starting from a sam- 
pling of possible nuclear fields A\, A\ or initial 
charge densities pf, pf in the colliding nuclei. Here 
we will present a calculation of the expectation 



4 



value of the flow fields carried out in the MV model 
with slowly varyin g av erage charge density fi. It 
has been shown in [19[ that the expectation value 
of the initial energy density is given by the prod- 
uct of the average source densities \xi in both 
nuclei, 
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ln 2 ^ Ml (r> 2 (f). (15) 



For the transverse Poynting vector let us discuss 
our expectations before calculating the result. First 
we note that the non-abelian terms from the co- 
variant derivatives would lead to expectation values 
of three-gluon correlators (AAA) in one of the nu- 
clei which under the assumptions of the McLerran- 
Venugopalan model vanish. Furthermore we can 
convince ourselves that because of the relation be- 
tween Eq and Bq given by A\ and A\ the two 
terms in Q are the same up to a sign and thus 
add up constructively. One can also easily argue 
that for dimensional reasons /iiV l /i2 + M2VV1 an d 
M1VV2 — A^V^ii and their contractions with e* J ' 
are the only four transverse vectors available after 
averaging. 

The first option, symmetric in the two nuclei and 
circulation-free, is picked by the 77-even flow term 



S+ = --cosh v — — 



In 2 -^-V 4 (fJ-ifJ-2) 



--coshr/V l e . (16) 



It mimics hydrodynamic flow. For the ?7-odd flow 
term we need to know the expectation value of two 
gluon fields in light cone gauge with one derivative. 
Following the same recipe that leads to the expres- 
sion d2|) for the gluon distribution and using our 
smoothness condition on /x(r) we obtain 



(( 



<1«Ai) (S x )A{(x ± )) = — 
Q 2 



(N? - 1 



-Sab 



x In ^-S7 l fj,(x±) (5 jl 6 ik - S il S jk - S kl 5 lJ ) (17) 

as the leading term contributing to S' l _. We refer 
the reader to a detailed calculation of this correla- 
tion function in a forthcoming publication [2lj |. It 
is then straight forward to obtain the rapidity-odd 
flow as 

ni r g 6 N c 



x In 
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Figure 4: Flow field V (black arrows) and energy density 
eo (shading) in the transverse plane for Au+Au collisions at 
6 = 6 fm. The nucleus centered at x = 3 fm travels into the 
plane which is the positive rj-direction. Left Panel: r) = 0. 
Right Panel: rj = 1. 




Figure 5: Same as Fig. [4] but plotted in the r\ — x-plane 
defined by y = 0. The flow will lead to a tilted fireball. 



5. Discussion of Results 

To summarize, the prediction of classical QCD 
for the initial average transverse energy flow nor- 
malized by the average initial energy density is 



V 1 = 



-- fcosh^^ 



- sinh 77 



/i 2 V>l - ^lV>2 

M1M2 



(19) 



which is independent of the UV cutoff and the IR 
regulator. In the following we have calculated V 1 
in several situations using Woods-Saxon profiles for 
incident nuclei. Fig. [4] shows the average flow field 
V 1 for the collision of two gold nuclei at impact 
parameter b = 6 fm in the transverse plane for two 
space-time rapidities. One can clearly see the evolu- 
tion in rapidity from a hydro-like flow field at r\ = 
to a preferred flow direction at forward rapidity. 



5 




Figure 6: The same as Fig. [5] for Au+Cu (Au traveling to 
the right). Left Panel: 6 = fm. Right Panel: 6 = 2 fm. 
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Figure 7: The same as Fig. [4] for Au+Cu at 6 = 2 fm. Left 
Panel: r) = 1, Right Panel: r] = —1. 



In Fig.[5]the same collision is shown in the r\ — x- 
plane. Clearly the flow tilts the fireball clockwise. 
The orientation of rotation is as if the gluon flux 
tubes preferred to expand in the wake of spectator 
nucleons in such a way that the flow increases with 
increasing separation from the spectators in rapid- 
ity. However this can not be taken literally as the 
origin of the effect. Our calculation is based on a 
small time expansion and the response of the en- 
ergy density to the flow will come in at the next 
order. Note that the normalizations of the vector 
fields in the figures are arbitrary. Typical values of 
the flow V 1 at the surface for Au+Au collisions is 
~ 0.1 at t = 0.1 fm/c at midrapidity. Fig. (Bjshows 
the average flow fields V 1 for Au+Cu collisions at 
impact parameters b = fm and b = 2 fm in the 
r\ — x-plane. In the central case the flow field leads 
to an expansion which is much more pronounced 
on the Cu-side of the system, consistent with the 
rule of thumb that flux tubes like to expand into 
the wake of spectators (which here are solely from 
the gold nucleus). The flow pattern becomes more 



involved for Au+Cu collisions at finite impact pa- 
rameter. In Fig. [7] the flow in the transverse plane 
is shown for forward and backward rapidity for the 
b = 2 fm Au+Cu system. We notice that the az- 
imuthal modulation of the flow is non-trivial but 
can again be understood through the position of 
spectator nucleons from the Au nucleus (centered 
at x = 1 fm) . This and the previous figures make it 
clear that 51. contributes not only to directed flow 
but also to the elliptic flow. 

The flow of energy in the classical field before Tth 
will translate into a flow of energy in the hydrody- 
namic phase after thermalization due to local en- 
ergy and momentum conservation [7|. One expects 
remnants of this flow to survive in hydrodynamics 
due to the inertia of fluid cells. In particular, this 
should result in a directed flow of particles which is 
odd in momentum rapidity y. In fact such a y-odd 
directed flow, measured by the first Fourier compo- 
nent vi , has been observed at RHIC 23|, \2M ■ The 
sign of the effect is consistent with the expectation 
from color glass, moreover the data points as a func- 
tion of rapidity could be fitted with a sinh y-shaped 
function. At this point it is too early to draw strong 
conclusions but the coincidence of sign and shape 
of the effect with data is encouraging. 

Some of the qualitative features of the flow field 
discussed here have been generated in hydrody- 
namic simulations by initializing a tilted source [25| , 

Our 
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e.g. postulated in the fire streak model 
calculation suggests that color glass could account 
for this phenomenon without invoking additional 
model assumptions. In addition, classical QCD 
adds several unique predictions in particular for the 
case of collisions of asymmetric nuclei. A system- 
atic study of flow as a function of rapidity and dif- 
ferent nuclear systems could find this unique finger- 
print of color glass. 
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